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ABSTRACT 

We study an accretion flow during the gravitational-wave driven evolution of 
binary massive black holes. After the binary orbit decays due to an interaction 
with a massive circumbinary disk, the binary is decoupled from the circumbinary 
disk because the orbital-decay timescale due to emission of gravitational wave 
becomes shorter than the viscous timescale evaluated at the inner edge of cir- 
cumbinary disk. During the subsequent evolution, the accretion disk, which is 
truncated at the tidal radius because of the tidal torque, also shrinks as the or- 
bital decay. Assuming that the disk mass changed by this process is all accreted, 
the disk becomes radiatively inefficient when the semi-major axis is several hun- 
dred Schwarzschild radii. The high-energy radiations, in spite of a low bolometric 
luminosity, are emitted from an accretion disk around each black hole long be- 
fore the black hole coalescence as well as the gravitational wave signals. The 
synchrotron process can notably produce potentially observable radio emissions 
at large distances if there is a strong, dipole magnetic field around each black 
hole. In unequal mass-ratio binaries, step-like light variations are seen in the 
observed light curve because the luminosity is higher and its duration time are 
shorter in the radio emission by the disk around the secondary black hole than 
those of the primary black hole. Such a precursor would be unique to not a single 
black hole system but a binary black hole system, and implies that binary black 
holes finally merge without accretion disks. 



Subject headings: black hole physics - accretion, accretion disks - gravitational 
waves - galaxies: evolution - galaxies: active - galaxies: nuclei - quasars: general 
- binaries: general 
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Introduction 



Astrophysical disks are ubiquitous and sources of active phenomena in the various sys- 
tem of the Universe; the star-compact object systems, star-planet systems, active galactic nu- 
clei ( AGNs), and so forth. Although the standard theory of accretion disks flShakura Sz Sunyaev 
19731 ) has succeeded in explaining them, there are still many unsolved problems. 



High-energy emission from galactic black hole binaries and AGNs cannot be produced by 
the accretion flow, except for the accretion flow accompanied by a corona, derived from the 
standard disk theory, where radiative cooling is so efficient that the thermal energy heated 
up by the viscosity is locally radiated away. If the radiative cooling becomes inefficient, 
the thermal energy is transported inward with accretion, and then the flow becomes hot. 
Such a flow can produce high-ener gy emission and is generally called as radiatively inefficient 
accretion flow (RIAF) (Chap. 9 of kato et allboosf ). 



The proto-ty pe model of the RIAF is an opt i cally thin advection-dominate d accre- 



tion flow ( ADAF) (llchimaru 



Chap. 9 of iKato et al. 



2008 



19771 : lNaravanfcYilll994 119951 : lAbramowicz et al.lll995t see also 
for a review). A low luminosity with high-energy emission can 
be explained by the ADAF with a significantly lower accretion rate than Me = Le/c 2 ~ 
2x lO^M^Moyr^ 1 ], where L E = 4vrcG'M bh / K es ~ 1.3 x lO^il^ergs- 1 ], c, G, M 7 , and k cs are 
the Eddington luminosity, light speed, gravitational constant, black hole mass normalized by 
1O 7 M , and electron scattering opacity, resp ectively. Observed pow er-low spectra of X-rays 
from Sgr A* was well reproduced by ADAFsf lManmoto et al.lll997l ). 



Most galaxies are thought to have massive black holes at their centers ( Kormendv Sz Richstone 



19951) and to coevolve with them( Magorrian et al.fl998 ; FerraresefcMerritt 



2000 



Gebhardt et al 



200of l. Galaxy merger leads to the mass inflow to the central region and then a nucleus of 
the merged galaxy is activated and black hole grows by gas accretion. During a sequence 
of processes, binary massive black holes with a sub-parsec scale separatio n are inevitably 



form e d before two black holes merge by emitting gravitation al radiation tjBegelman et al. 



1980; Escala et al 



gaseous disks (llvanov et al 



2009 



2005; 



Dotti et al. 



2007 



rounding star s (IMerritt Sz Milosavlievic 



Cuadra et al. 



2009 



Mayer et al.l 120071 ) . B y interaction with sur 



20051 and references therein; Berentzen et al 



Haiman et al. 



20091). other m assive black holes (llwasawa et al. 



2009) 



1999l:lGould fc Rix |[200o l; lArmitage Sz Nataraian II2OO2I; lHava"saki 



20061 ). and infalling dwarf galaxiesf lMatsui Sz Habd 120091 ) . the binary evolves towards co- 
alescence within Hubble time, although whether the binary finally coalesces still remains 
uncertainty theoretically (e.g.. lLodato et al.l 120091 ) and observationally. Then, gravitational 
radiation emitted prior to coalescence is detectable with the Laser Interferometer Space 
Antenna (LISA) and Pulser Timing Arrays(PTAs). 
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An electromagnetic counterpart of gravitational wave signal plays a key role in determin- 
ing the redshift of the source. By combining the redshift with its luminosity distance the grav- 
itational wave determines, one can obt ain the distance-redshift relation, which is one of best 
observational probes of the dark energy fjSchutzlll986l ; lHolzfeHughesll2005l) . Some types of the 



electromagnetic counte r parts have been studied such as afterglows ( iMilosavlj evicfcPhinney 



2005; Rossi et al. 2010 



prompt e missionsdChang et al. 



Corrales et al. 
200_9T 



emissions (IHavasaki et al 



200 



20081 : IHavasaki fc Okazaki 1 120091 ) . 



7 



2010 



Bode et al. 



TanakafeMenoul 120101) and pre cursors as the 



201p[ Istone fe Loebl l2010h and periodic 



20081 ; iBogdanovic et al. 



2008 



MacFadyen fc Milosavljevic 



In this Letter, we investigate RIAFs triggered by the rapid orbital decay due to emis- 
sion of gravitational wave. The letter is organized as follows. The evolution of binary 
massive black holes interacting with a massive circumbinary disk are briefly described in 
Section 2. In Section 3, we derive the possible accretion flows and their luminosities during 
the gravitational-wave driven evolution. Section 4 is devoted to summary and discussion. 



2. Final Parsec Evolution of Binary Massive Black Holes 



We first briefly describe the evolution of binary massive black holes interacting with a 
massive circumbinary disk in the framework of coevol ution of massive black holes and their 
host galaxies. The detailed description can be seen in IHavasaki et al.l ( 120101 ). 



Binary black holes are considered mainly to evolve via three stages (IBegelman et al. 



19801 ). Firstly, each of black holes sinks independently towards the center of the common 
gravitational potential due to the dynamical friction with neighboring stars. When the 
separation between two black holes becomes less than 1 pc or so, angular momentum loss 
by the dynamical friction slows down due to the loss-cone effect and a massive hard bi- 
nary is formed. This is the second stage. The binary hardens at the radius where the 
kinetic energy per unit mass of the star with the velocity d ispersion equals to the bind- 
ing energy per unit mass of binary black holes ( iQuinlanl Il996f ) . Its hardening radius is de- 
fined as ah ~ 8.5 x 10~ 1 [pc](q/(l + q) 2 )M* 2/Pa , where q = M s /M p is the ratio of primary 
black-hole mass, M p , to secondary black-hole mass, M s . Here, it is assumed that the total 
black-hole mass is tightly correlated with the velocity dispersion: M 7 = f3\ (cr ;iI /200kms~ 1 ) /32 
with B-\ = 16.6 and 02 = 4.86, whe re a* shows the one- dimensional velocity dispersion of 



stars fjMerritt fc Milosavljevic II2005I ). Finally, the semi- major axis of the binary decreases 
to the radius at which the gravitational radiation dominates, and then a pair of black holes 
merge into a single, more massive one. 




Fig. 1. — (a)Black-hole mass dependence of characteristic radii of the accretion disk around 
the primary black hole in the binary with q = 0.1. The solid line represents the tidal radius 
when the binary is decoupled from the circumbinary disk. The dashed line shows the critical 
transition radius, r tjC , where the disk-outer edge corresponds to the transition radius, r t , 
from the standard disk to the RIAF. The dotted line shows the radius, rb, , where the disk- 
outer edge corresponds to the boundary, r b , between the inner, radiation-pressure dominated 
region and the middle, gas-pressure dominated region. (b)Evolution of characteristic radii 
for Mbh = 1O 7 M . The solid line, dashed line, and dotted line represent the disk-outer edge, 
r Q (a), r t , and rb, respectively. The intersection between the solid line and the dashed line 
shows r t)C , whereas the one between the solid line and the dotted line shows r^o- 

2.1. Triple-Disk Evolution 

The circumbinary disk would be formed around hardening of the binary. The inner 
edge of circumbinary disk, r; n , is then defined « 2a ~ 1.8[pc](g/(l + g) 2 )M 7 1 ~ 2//32 (a/a h ) 
, where a is the semi-major axis of binary. For simplicity, the circumbinary disk is as- 
sumed to be a steady, axisymmetric, geometrically thin, self-regulated, self-gravitating but 
non-fragmenting with a Keplerian rotation and accretion rate, M acc = [r]/ e)M^di with the 
Eddington ratio, ?7, mass-to-energy conversion efficiency, e, and Shakura-Sunyaev type vis- 
cosity parameter, a sg . 

The circumbinary disk and binary exchanges the energy and angular momentum through 
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Fig. 2. — Black-hole mass dependence of the normalized accretion rate and corresponding 
normalized luminosity in the binary with q = 0.1. The normalized accretion rate, ttiriaf j 
is evaluated at the critical transition radius, r tx . The solid line and dashed line show ttiriaf 
and corresponding normalized luminosity, L mj ^ F /L E , respectively. The dotted line shows the 
normalized critical accretion rate, fn c . 



the tidal/resonant interaction, which leads to the orbital decay of the binary. On the 

other hand, the gas at th e inner edge of circumbinary disk overflows onto the central 

binar v llHayasaki et al 1 l2007h . Then, an accretion disk is formed around each black hole (IHayasaki et al. 



20081 ) . The mass tran sfer therefore adds its angular momentum to the binary via accretion 
disks (IHayasaki II2009I ). In a steady state, the mass transfer rate equals t o M arr . Since it i s 
much smaller than the critical transfer rate defined by equation (41) of IHayasaki I (120091 ). 
the additional torque to the binary by the mass transfer can be neglected. The orbital- 
decay timescale is then given as t gas ~ 3.1 x 10 8 [yr](g/(l + <?) 2 )?7aieo.i, where 770.1 = 77/0.1 
and eo.i = e/0.1. The total mass supplied to the central binary during the orbital decay is 



obtained as 



Q 



M sup = M acc t gas ~ 6.2 x lO b [M ]M 7 ^— 



M sup thus gives an upper limit of mass of two accretion disks. 



2.1.1. Decoupling of Circumbinary Disk 

The orbital -decay timescale of the binary by gravitational-wave emission can be written 
by iPeters! (119641 ) as 



t 



?) 2 (1 



3 2\7/2 



q 



/(e) 



(2) 



where e is the orbital eccentricity, /(e) = (1 + 73e 2 /24 + 37e 4 /96) is a function of the orbital 
eccentricity, and r g = 2GM\^/c 2 is the Schwarzschild radius of total black hole mass. For 
simplicity, the orbital eccentricity set to be zero, unless otherwise noted, in what follows. 

As the binary evolves, the circumbinary disk behaves as the standard disk in the 
inner region for < 10 7 M Q , whereas the disk still remains to be self-gravitating for 

> 10 7 M Q . When the orbital-decay timescale is shorter than the viscous timescale mea- 
sured at the inner edge of the circumbinary disk, the binary is decoupled from the circumbi- 
nary disk. For the self-gravitating disk, the decoupling radius can be defined as a&/r g = 
[(8/5)(ct vis Jr s )(q/(l + q) 2 )] 2/7 , where ~ 5.6 x ^b^ g/O-Oe)- 1 / 3 ^ eg^M^ 3 is 

the viscous timescale measured at a = r g (lHayasaki et al.ll2010l ). For the standard disk, 



l/4 £ l/4 M 5/4 
-0.1 m 7 



a d /r g = [(8/5)(ct viSig /r g )(g/(l + g) 2 )] 4 / 11 , where t vis , g ~ 5.9 x 10 2 [yr](a ss /0.1)- 1 r / o 
After the binary is decoupled from the circumbinary disk, the system consists of only an ac- 
cretion disk around each black hole. 



3. Subsequent Double-Disk Evolution 

For simplicity, we focus on the accretion disk around the primary black hole in what 
follows. Assuming that the binary has a circular orbit, the disk outer edge can be defined 
by the tidal radius as r G (a) ~ 0.9^ ~ 4.2 x KT 1 (M p /M bh ) 1/3 a, where M hh = M v + M f 



and i?i shows the inner Roche radius of the primary black hole for q > O.lfe.g.. iFrank et al. 



20021 ). The logarithmic differentiation of the tidal radius is given by 
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where M p and M s are the growth rate of primary black hole and that of secondary black 
hole, respectively. M p /M p <C 1 and M s /M s <C 1 are held during the evolution. 

The accretion disk, which is truncated at the tidal radius because of its tidal torque, 
shrinks with the orbital decay due to gravitational-wave emission. The disk mass lying 
between r + Ar and Ar can be then estimated as 

/•r+Ar 

AM = 2vr / S(r)rrfr « 27rrS(r)Ar. (4) 



We assume that the disk mass changed by shrinking of the tidal radius steadily accretes 
onto the black hole without being emitted as jet or wind. Combining equation (j3J) with 
equations (J2J) and Q, we obtain the accretion rate normalized by M E : 

2vrr 2 £r 2vrr 2 S(r) 1 
m = —. = 5 

M E r tgw M E 



Following IShakura fc Sunyaevl (119731 ). the standard disk structure has three distinct 
regions for a given accretion rate due to the source of opacity and pressure. In the inner 
region, the pressure and opacity of the disk are dominated by the radiation pressure and 
electron scattering, respectively. The pressure and opacity are dominated by the gas pressure 
and opacity of the electron scattering, respectively, in the middle region, whereas they are 
dominated by the gas pressure and opacity of the free-free absorption, respectively, in the 
outer region. At the present case, the outer region is neglected because the disk-outer edge is 
always less than the boundary between the middle region a nd the oute r regio n. The surface 



density in the inner region and middle region is written by iKato et al.l (120081 ) as 



1 1 \ S mid (r/r g )- 3 / 5 m 3 / 5 r > r b ' 1 } 

where E in = 1.0xl0 3 a a J(l+g) 3 / 2 /^ 1 fe/cm 2 ], S mid = 4.3x 10 31 / 5 a~i /5 M 7 1/5 (l+g)- 4 / 5 / 3/5 [g/cm 2 ] 
with / — 1 — A/3r g /r, and 



n, 



1.8 x 10 n / 7 KiM 7 ) 2 / 21 m 16 / 21 (l + g)" 23 / 21 (7) 



is the boundary between the inner region and the middle region. 

A mass of accretion disk around the primary black hole at the decoupling is defined by 
M a d = S(r)r<ir, where r isco = 3r K /(l + o) shows the radius of innermost stable circular 

J f*isco o/ \ / 

orbit around the Schwarzschild black hole. Since M a d is an increasing function of black hole 
mass and M ad /M bh < 10~ 3 for M bh = 1O 8 M and q — 1, M ad is always less than M sup over 
the whole mass range. 
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From equations (El) and ([6]), we finally obtain the normalized accretion rate as 



m 



(27rHE in /t gw )M E (r (a)/r g ) 



1/2 



v7/4 



r < r b , 



5/2 



(27rr 2 £ mid /t gw )/M E (r G (a)/r g ) 7 / 2 r > r b>0 



where r^o is a newly defined boundary between the inner region and middle region, which 
can be written as 



rv 



16tt 



Q 



r„c 



5 (l + g) 2 M E 



8/25 y2/5 
mid 



,2/25 ' 



(9) 



Note that the entire region of the disk is radiation-pressure dominated when r (a) < r^o- 

The viscous timescale of the accretion disk around the primary black hole is written by 
t vis pa y^GM hh r/(l + g)/(a;g S c 2 ), where c s is the sound velocity of the middle region of the 
disk. The ratio of t vis to t gw is then given by 



^~(l + g) 

fcgw 



7/5 

r (a),' 



(10) 



where equation ([8]) is substituted in the derivation process. Note that this ratio is indepen- 
dent of the viscosity parameter and black hole mass. The assumption that the disk is the 
steady state during the gravitational- wave driven evolution is thus justified in most of region 
of the disk because of t v i S < t gw . From equation f lTUj) . t vis is slightly longer than t gw at the 
disk-outer edge. If the outer region of the disk is, however, hotter by the tidal dissipation, 
t v i S would become shorter than t gw . In order to confirm whether this solution is realized, the 
time-dependent behavior including the effect of the tidal dissipation will be investigated in 
a subsequent paper. 

Next, we consider that the disk becomes a bimodal system of standard disk-RIAF, where 
the standard disk in the cool outer region transits to the RIAF in the inner hot re gion. We 
emplo y a transition radius from the standard disk to the RIAF, which is defined by 
jl996h as 



Homma 



-2 4 -2 • 

^o.i e o.i r " 



(ii) 



The critical transition radius where the disk-outer edge corresponds to the transition radius 
is then written by using equation (j2J) and (JSJ) as 



^ ~ 9.0 x 10" 



14 -2 
a 0.1 e 0.1" 



q 



1 + g) 43 / 18 



r g cS m i(j 
Mr 



5/12 



(12) 
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Note that the disk becomes the bimodal system when r D (a) < r t> , 



Fig. [T] shows the mass dependence of characteristic radii of the accretion disk in the 
binary with q = 0.1, and their evolution for = 1O 7 M . All the radii are normalized 
by r g . In panel (a), the solid line is the tidal radius when the binary is decoupled from the 
circumbinary disk. The dashed line and dotted line are r tjC and rb, OJ respectively. Panel (a) 
shows rb j0 < J"t )C < ^o(od) in the range from 10 5 M Q to 1O 8 M . In panel (b), the solid line, 
dashed line, and dotted line are r D (a), r t , and r^, respectively. The intersection between the 
solid line and the dashed line represents r t)C , whereas the one between the solid line and the 
dotted line shows r^ Q . It is noted from panel (b) that the disk becomes the RIAF before the 
disk is radiation-pressure dominated. 

The l uminosity of the RIAF normalized by the Eddington value can be written as(cf. 
Chap. 9 of lKato et al.ll2008l ) L/L E « e (m/m c ) rh, where m c ~ 1.0 x lO^e^Ja 2 , A is the critical 
accretion rate normalized b y Me. Note that there are no possible solutions of RIAFs for 
m > rh c (|NarayanfeYilll995l ). Fig. [2] shows the mass dependence of the normalized accretion 
rate and corresponding normalized luminosity in the binary with q = 0.1. The solid line and 
dashed line show the normalized accretion rate evaluated at the critical transition radius, 
^riaf, and corresponding normalized luminosity, Lrjaf/^e, respectively. The dotted line 
shows rh c . It is clear from the figure that there is a possible solution for the RIAF because 
of rh c > ?Ti riaf- Both m M AF ~ 2.0 — 3.0 x 10~ 2 and L R iaf/£e ~ 5.0 — 9.0 x 10" 4 slightly 
decrease with the black hole mass in the range from 1O 5 M to 1O 8 M . 

The virial temperature evaluated at the critical transition radius, T v - ir ~ GMbhm p /r tiC , 
is 10 10 ~ 11 K in the current-mass range. Such a high temperature disk can produce the X- 
ray emissions by the inverse Compton scattering and thermal bremsstrahlung, and radio 
emission by the synchrotron emission process. The luminosity of the synchrotron emiss ion 
at the critical transition radius can be estimated by (Chap. 6 of lRybickifcLightmanlll979l ) 



^sync 



2 a T (3 2 



9 M E 1 - 



3.1 x 10 11/2 



%1«0.1 



3/2 M7 -3/2 (1+g) l/4 



Hi 



-17/4 



(13) 



where f3 2 = (r t)C /r g ) 1 (l+q) 1 , the electron number density is evaluated by n c ~ S(r tjC ) / (m p H), 
where H is the disk-scale hight, and mag netic field is evaluated by B ~ 2.7x 10 5 v /t] , 1 /M 7 (1 + 
?) 1 ^ 2 ( r t,c/ r g)~ 3 , based on an assumptio n of a dipole field wit h the strength of an equipartition 
value around a non-rotating black hole (jDemerfcFinkell2008l ). For the binary with 1O 7 M and 
q = 0.1, L sync /LE ~ 5.0 x 10~ 2 , which is significantly larger than that of other processes 
in the current case. This is a promising candidate as an electromagnetic counterpart of 
gravitational wave emitted at the black hole coalescence. 
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In an unequal mass-ratio binary, the critical transition radius of the disk around the 
primary black hole becomes larger than that of the secondary black hole. This difference 
causes step-like light variations in the radio emissions. Since the duration time approximately 
equals to the dynamical time evaluated at the critical transition radius, the difference of 
duration time can be written by Aeration = V / 2(^g/c)(r tjC /r g ) 3 / 2 (l + g) 1 ^ 2 (l — g 5 ^ 6 ). For the 
binary with Mbh = 1O 7 M and q = 0.1, the difference of duration time becomes Aeration ~ 
2.0 x 10 -3 yr. On the other hand, the synchrotron-luminosity ratio of the secondary black 
hole to that of the primary black hole is £ synCjS /L synCiP — q~ 5 ^ 3 . For the binary with q = 0.1, 
the emission from the secondary black hole shows one-order of magnitude higher luminosity 
but shorter duration time than those of the primary black hole. The total light curve would 
thus exhibit step-like variations. 



4. Summary&Discussion 

We study radiatively inefficient accretion flows onto each black hole during the gravitational- 
wave emission driven evolution of binary massive black holes in the framework of coevolution 
of massive black holes and their host galaxies. 

The binary system initially consists of the accretion disk around each black hole and 
the massive circumbinary disk surrounding them. After the binary is decoupled from the 
circumbinary disk, the outer edge of each accretion disk shrinks with the orbital decay due 
to emission of gravitational wave. Assuming that the disk mass changed by the rapid orbital 
decay is all accreted, the accretion flow becomes radiatively inefficient when the disk-outer 
edge is less than the critical transition radius. 

The bolometric luminosity corresponding to the accretion rate estimated at the critical 
radius is 5.0 - 9.0 x 10 _4 L E in the range from 10 5 M© to 10 8 M Q . The RIAF produces X- 
ray emissions by the inverse Compton scattering and thermal bremsstrahlung, and radio 
emissions by the synchrotron process. Assuming that the magnetic field is a dipole field 
around a non-rotating black hole and its strength is the equipartition value, The radio 
emission by the synchrotron process is much more luminous than those of other processes 
and exhibits step-like variations by a sum of light curve from each component of the binary. 
Such a precursor implies that each accretion disk are drained because of the high energy 
radiation emitted by each disk prior to coalescence. 

The peak frequency of the synchrotron emission is given by v c = 37 2 e-Bsina/(47rm e c), 
and 7, m e , and a are the Lorenz factor, electron mass, and pitch angle between the magnetic 
field and the electron velocity, respectively. The peak- frequency ratio of the radiation emitted 
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by the disk around the secondary black hole to that of the primary black hole can then be 
estimated by u CB /u cp ~ g -3 / 2 because the Lorenz factor and magnetic field strength depend 
on the mass ratio. The spectrum is therefore expected to have a double peak when one 
observe in a wide range of radio wavelengths. 

A step-like nature of the light curve and double peak nature of the spectrum are unlikely 
to be seen in emissions from an accretion flow around a single massive black hole, and thus 
gives an unique photon diagnosis of binary black hole coalescences at large distances. 
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